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Abstract 

Background: Few data on the relationship between temperature variability and childhood pneumonia are 
available. This study attempted to fill this knowledge gap. 

Methods: A quasi-Poisson generalized linear regression model combined with a distributed lag non-linear model 
was used to quantify the impacts of diurnal temperature range (DTR) and temperature change between two 
neighbouring days (TCN) on emergency department visits (EDVs) for childhood pneumonia in Brisbane, from 2001 
to 2010, after controlling for possible confounders. 

Results: An adverse impact of TCN on EDVs for childhood pneumonia was observed, and the magnitude of this 
impact increased from the first five years (2001-2005) to the second five years (2006-2010). Children aged 5-14 
years, female children and Indigenous children were particularly vulnerable to TCN impact. However, there was no 
significant association between DTR and EDVs for childhood pneumonia. 

Conclusions: As climate change progresses, the days with unstable weather pattern are likely to increase. 
Parents and caregivers of children should be aware of the high risk of pneumonia posed by big TCN and take 
precautionary measures to protect children, especially those with a history of respiratory diseases, from climate 
impacts. 
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Background 

Pneumonia is the top cause of mortality in children 
under five years [1]. It is estimated that in 2010, world- 
wide, there were 120 million episodes of pneumonia in 
children younger than five [1]. Pneumonia is highly 
preventable, hence it is particularly important to 
explore the risk factors which drive the incidence of 
pneumonia and further to prevent children from being 
exposed to these risk factors. 

Many nutritional, socioeconomic and environmental 
factors are involved in the occurrence of pneumonia 
[2-4]. As climate change proceeds, the possible impact 
of climate factors on pneumonia transmission has 
attracted increasing research attention [4,5]. Both high 
and low temperatures have been reported to be associ- 
ated with increased pneumonia incidence [6,7]. How- 
ever, the potential impact of temperature variability on 
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childhood pneumonia has not been researched yet, 
though big temperature changes may influence the 
function of respiratory system [8]. 

There are several ways to define temperature variabil- 
ity [9]. For example, the difference in daily maximum 
and minimum temperatures (i.e., diurnal temperature 
range (DTR)) [10], and the mean temperature differ- 
ence from one day to the next (i.e., temperature change 
between two neighbouring days (TCN)) [11,12]. Previ- 
ous studies have highlighted that big DTR or TCN may 
affect the respiratory system of human [10-12], espe- 
cially for children [13]. We hypothesized that great 
DTR or TCN might be associated with increase in 
childhood pneumonia cases, and we used the data on 
emergency department visits (EDVs) for childhood 
pneumonia in Brisbane from 2001 to 2010 to test our 
hypothesis. 
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Methods 

Data collection 

Data on EDVs from 1 st January 2001 to 31 st December 
2010 classified according to the International Classification 
of Diseases, 9 th version and 10 th version (ICD 9 and 10) 
were supplied by Queensland Health. We extracted those 
cases coded as pneumonia (ICD 9 codes: 480-486; ICD 
10 codes: J12-J18) in children aged 0-14 years. Data on 
climate variables, including maximum and minimum tem- 
peratures, rainfall and relative humidity, were obtained 
from Australian Bureau of Meteorology. DTR was calcu- 
lated as daily maximum temperature minus daily minimum 
temperature [10]. Daily mean temperature was the average 
of daily maximum and minimum temperatures, and TCN 
was calculated as mean temperature of the current day 
minus mean temperature of the previous day [11]. Data 
on air pollutants, including daily average particular 
matter < 10 um (PM 10 ) (ug/m 3 ), daily average nitrogen 
dioxide (N0 2 ) (|ig/m 3 ) and daily average ozone (0 3 ) 
(ppb), were retrieved from the Queensland Department of 
Environment and Heritage Protection. Ethical approval was 
obtained from the Human Research Ethics Committee of 
Queensland University of Technology (Australia) prior to 
the data collection (number: 1000001168). 

Data analysis 

Distributed lag non-linear model (DLNM) was devel- 
oped to incorporate both lagged and the non-linear 



effects of temperature on mortality or morbidity [14,15]. 
Previous studies have revealed that there might be a 
lagged effect of temperature variability on human health, 
and the relationship between temperature variability and 
respiratory diseases appears to be non-linear [11-13,15]. 
Thus, we used DLNM to incorporate the non-linear and 
lagged effect [14]. A quasi-Poisson generalized linear 
regression combined with DLNM was used to quantify 
the association between DTR (or TCN) and EDVs for 
childhood pneumonia. 
DTR model: 

Y t ~ quasiPoisson(fi t ) 

Log (ji t ) = a + pDTR tj i + ns{T th 3) + ns(RH u 3) 

+ns(PM m , 3) + ns(0 3t , 3) + ns{N0 2tl 3) 
+ns(Time tl 8) + v\ x Holiday + rj 2 DOW t 

TCN model: 

Y t ~ quasiPoisson(fi t ) 

Log (p t ) = a+pTCNtji + ns{T th 3) + ns(RH t , 3) 

+ns(PM 10t , 3) + ns{0 3tl 3) + ns{N0 2t , 3) 
+ns(Time tl 8) + r\ x Holiday + rj 2 DOW t 

Where t is the day of the observation; Y t is the 
observed daily childhood pneumonia on day t; a is 
the model intercept; DTR t j is a matrix obtained by 
applying DLNM to DTR, TCN t j is a matrix 
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Table 1 Spearman's correlation between daily weather conditions, air pollutants and pediatric pneumonia in Brisbane, 
Australia, from 2001-2010 





Mean temperature 
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Mean temperature 
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-0.44* 


1.00 
















TCN + 


0.14* 


0.03 


1.00 














Relative humidity 


-0.06* 


-0.33* 


0.09* 


1.00 












Rainfall 


0.16* 


-0.49* 


-0.03 


0.38* 


1.00 










PM 10 


0.16* 


0.31* 


0.08* 


-0.29* 


-0.31* 


1.00 








o 3 


0.03 


0.20* 


0.03 


-0.28* 


-0.07* 


0.31* 


1.00 






N0 2 


-0.67* 


0.42* 


-0.02 


0.13* 


-0.14* 


0.01 


-0.10* 


1.00 




Pneumonia 


-0.37* 


0.21* 


-0.01 


0.03 


0.06* 


0.07* 


0.10* 


0.28* 


1.00 



*P < 0.01; + DTR, Diurnal temperature range; TCN, temperature change between two neighbouring days. 



obtained by applying DLNM to TCN, T t>i is a matrix ob- 
tained by applying the DLNM to temperature; ft is vector of 
coefficients for T t>i , and / is the lag days; ns(RH p 3) is a nat- 
ural cubic spline with three degrees of freedom for relative 
humidity; ns(PM 10p 3) is a natural cubic spline with three 
degrees of freedom for PM 10 ; ns(0 3p 3) is a natural cubic 
spline with three degrees of freedom for 0 3 ; ns(N0 2p 3) is a 
natural cubic spline with three degrees of freedom for N0 2 ; 
ns(Time p 8) is a natural cubic spline with eight degrees of 
freedom per year for long-term trend and seasonality; 
Holiday is the public holiday, and DOW t is the categor- 
ical day of the week with a reference day of Sunday. 

Specifically, DTR (or TCN) and lag were incorporated 
using a "natural cubic spline-natural cubic spline" 
approach. The model included lags up to 21 days for 
DTR (or TCN) and mean temperature [16]. We used 
lags up to 10 days, for all other confounders (i.e., relative 
humidity, PM 10 , 0 3 , and N0 2 ). 

All data analysis was conducted using the R statistical 
environment (v 2.15), and "dlnm" package was used to 
fit the regression. In the sensitivity analysis, we changed 



the df for DTR, TCN and time. We also excluded 2009 
data as there was a big pneumonia spike in 2009. 

Results 

Figure 1 shows the time-series distributions of childhood 
pneumonia, mean temperature, DTR and TCN, revealing 
that there was a seasonal trend in childhood pneumonia, 
mean temperature, and DTR. The 2009 pneumonia peak 
(due to H1N1 flu pandemic) is also revealed in this 
figure. To explore the crude relationship between each 
variable, we calculated the Spearman correlations between 
climate variables, air pollutants and childhood pneumonia 
(Table 1). There was a negative correlation between DTR 
and mean temperature (r = -0.44, P < 0.01). TCN was 
positively correlated with mean temperature (r = 0.14, 
P < 0.01). Further, DTR was positively correlated with 
childhood pneumonia (r = 0.21, P < 0.01), while no signifi- 
cant correlation between TCN and childhood pneumonia 
was observed. No correlation coefficient was greater than 
0.5, meaning that multi-collinearity is unlikely a big issue 
in the subsequent modelling. 




i 1 1 1 1 1 1 i i i i 

0 5 10 15 20 25 30 -10 -5 0 5 



DTR (°C) TCN (°C) 

Figure 2 The overall effects of DTR and TCN on pediatric pneumonia in Brisbane, from 2001 to 2010. 
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Figure 2 shows the exposure-response relationship 
between temperature variability and childhood pneu- 
monia (the modelling results). No significant associ- 
ation between DTR and childhood pneumonia was 
observed. In contrast, a big temperature decrease from 
one day to the next (TCN < -2°C), increased the risk of 
childhood pneumonia. Since TCN < -2°C is associated 
with an increase in childhood pneumonia, we subse- 
quently calculated the number of days with TCN < -2°C 
every year. There were more than 50 days with TCN 
below -2°C every year, with most days with temperature 
drop > 2°C occurring in the second half of each year (June 
to December) (Figure 3). Figure 4 shows the pattern of 
lagged effects of TCN on childhood pneumonia, revealing 
that TCN effect lasted for nearly three weeks. Figure 5 
depicts that older children (5-14 years vs. <5 years), 
female children (vs. male), and Indigenous children 
(vs. non-indigenous) appeared to be more vulnerable 
to the TCN impact. 

As there was a distinct seasonality in childhood pneu- 
monia, with the peak in winter, we specifically analysed 
the TCN impact on childhood pneumonia in summer 
(December, January, and February) and winter (June, July 
and August), and found that this impact mainly occurred 
in winter (Figure 6). However, in summer an increased 
relative risk of childhood pneumonia was also detected 
as TCN >0°C, but it was not statistically significant. 

To test whether there was a change over time in the 
effect of TCN on childhood pneumonia, we splitted the 
ten years into two periods (2001-2005 and 2006-2010). 
Figure 7 reveals that the effect of TCN on childhood 
pneumonia during the 2 nd period was much greater than 
it was during the 1 st period. 

Figure 8 shows that after excluding the 2009 data, the 
magnitude of TCN effect on childhood pneumonia in 
Brisbane reduced, although the shape of the TCN- 



TCN effect 

o 

OJ _ . 
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Figure 4 The lagged effect of TCN on childhood pneumonia in 
Brisbane, from 2001 to 2010. 

V. J 

pneumonia relationship was similar. Figure 9 shows 
the effects of TCN on childhood pneumonia in different 
subgroups after excluding the 2009 data, revealing that 
subgroups vulnerable to TCN effect remained largely 
unchanged. We also compared the effect of TCN on 
childhood pneumonia between 2001-2005 and 2006- 
2010 (without 2009), and found that the TCN effect on 
childhood pneumonia during the 2 nd period was still 
greater than it was during the 1 st period (Figure 10). 

Discussion 

This study quantified the impacts of both DTR and 
TCN on childhood pneumonia and yielded several 
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Figure 5 The effect of TCN on the total-, age-, gender- and ethnic-specific childhood pneumonia in Brisbane, from 2001 to 2010. 



notable findings. A big temperature decrease from one 
day to the next (TCN < -2°C) may increase the EDVs 
for childhood pneumonia, and this effect lasted for 
around three weeks. Every year, there were more than 
50 days with big TCNs, and these big TCNs mainly oc- 
curred in winter. Children aged 5-14 years, female 
children and Indigenous children were particularly at 
risk. Further, there was a change in the effect of TCN 
on childhood pneumonia over time. No significant re- 
lationship between DTR and childhood pneumonia 
was observed. 

Children are particularly vulnerable to both extreme 
temperatures [17] and temperature variation [13], due 
partially to their relatively less-developed thermoregu- 
lation capability [18]. In this study, we found that a 
sharp temperature drop was followed by significantly 
increased EDVs for childhood pneumonia, and the 
TCN impact lasted for roughly three weeks. The time 



lag of TCN impact is probably due to two factors: (1) 
the temperature change which exceeds certain limits 
may take a few days to trigger subsequent symptoms in 
children with underlying conditions; (2) there may be 
another delay between the onset of symptoms and seek- 
ing for medical attention. 

Some studies have observed significantly increased 
respiratory-related mortality associated with big TCN 
[12] while others did not find significant effects [11]. 
Our study stands out of previous studies by specifically 
focusing on childhood pneumonia and controlling for 
a range of possible confounders. In this study, we also 
found age, gender and Indigenous-status modified the 
relationship between TCN and pneumonia. The school- 
aged children (5-14 years) were more sensitive to TCN 
compared with younger, which might be because they play 
outdoors more often and thus exposed more to the out- 
door temperature change. The difference in vulnerability 
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Figure 7 The overall effects of TCN on childhood pneumonia 
during two periods. 



to TCN between two genders may be due to their body 
composition [19], though some researchers argued that 
such an effect is variable among locations and populations 
[20]. Our results also suggest that Indigenous children 
were more sensitive to TCN effect compared with non- 
Indigenous children. Previous studies have reported 
that the burden of pneumonia in Indigenous children 
is 10 to 20 fold higher than non-Indigenous children, 
and they have longer hospital admissions and are more 
likely to have multiple admissions with pneumonia [21]. 
Most Indigenous children have limited access to infra- 
structures, and experience more poverty than non- 
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Figure 8 The overall effect of TCN on childhood pneumonia in 

Brisbane, from 2001 to 2010 (excluding 2009). 
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Indigenous children, possibly resulting in their greater 
vulnerability to TCN impact [22]. 

As climate change progresses, not only the global sur- 
face average temperature, but also the frequency of un- 
stable weather patterns (e.g., sharp increase/decrease in 
temperature) will increase [23], which poses a significant 
challenge to public health sectors. We found that the ef- 
fect of TCN on childhood pneumonia during 2006-2010 
was greater than it was during 2001-2005. This finding 
indicates that children might be vulnerable to sharp 
temperature decrease in the future if unstable weather 
patterns occur as projected. Elucidating the impact of 
temperature variability on children's health is essential 
for the improvement of public health. The findings of 
our study not only remind parents and children's care- 
givers to take good care of children in the days of big 
TCN, but also imply that government should take 
temperature variability into account while developing 
early warning systems for controlling and preventing 
childhood pneumonia. 

This study has two major strengths. First, this is, to 
our knowledge, the first study to look at the impact of 
DTR and TCN on childhood pneumonia. Second, the 
change over time in the TCN effect on childhood pneu- 
monia which we observed in this study may encourage 
future studies to explore the temporal variability of TCN 
impacts on children's health. Two weaknesses should be 
acknowledged. First, this is a one city study, which means 
further the interpretation of our findings should be cau- 
tious. Second, to some extent, biases in exposure and/or 
outcome measures may be inevitable because we used the 
aggregated data on temperature and EDVs for childhood 
pneumonia. 
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Figure 9 The effect of TCN on the total-, age-, gender- and ethnic-specific childhood pneumonia in Brisbane, from 2001 to 2010 
(excluding 2009). 
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Figure 10 The overall effects of TCN on childhood pneumonia 
during two periods (excluding 2009). 



Conclusions 

No relationship between DTR and childhood pneumonia 
was observed. A sharp temperature decrease from one 
day to the next had an adverse impact on childhood 
pneumonia, especially in winter, and the magnitude of 
this impact increased in recent years. In addition, that 
effect of TCN on the risk of contracting pneumonia ap- 
peared to differ by age, gender, and race/ethnicity. As 
climate change continues, unstable weather patterns 
may be more frequent, and the findings of this study 
have important implications for public health policy to 
protect children from the impact of pneumonia. 
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